GLOBAL CLASSICAL SOLUTIONS OF THE RELATIVISTIC 

VLASOV-DARWIN SYSTEM WITH SMALL CAUCHY DATA: 

THE GENERALIZED VARIABLES APPROACH 

REINEL SOSPEDRA-ALFONSO, MARTIAL AGUEH, AND REINHARD ILLNER 
O 

' Abstract. We show that a smooth, small enough Cauchy datum launches 

a unique classical solution of the relativistic Vlasov-Darwin (RVD) system 
globally in time. A similar result is claimed in [15| following the work in |13| . 
Our proof does not require estimates derived from the conservation of the total 
l/i ' energy, nor those previously given on the transverse component of the electric 

field. These estimates are crucial in the references cited above. Instead, we 
exploit the formulation of the RVD system in terms of the generalized space 
and momentum variables. By doing so, we produce a simple a-priori estimate 
on the transverse component of the electric field. We widen the functional 
(~H , space required for the Cauchy datum to extend the solution globally in time, 

and we improve decay estimates given in |15| on the electromagnetic field and 
its space derivatives. Our method extends the constructive proof presented in 
|14) to solve the Cauchy problem for the Vlasov-Poisson system with a small 
initial datum. 
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The relativistic Vlasov-Darwin (RVD) system can be obtained from the Vlasov- 
Maxwell system by neglecting the transverse component of the displacement current 
in the Maxwell- Ampere equation. Precisely, consider an ensemble of single species 
charged particles interacting through the self-induced electromagnetic field. Let 
f{t,x,p) denote the number of particles per unit volume of the phase-space at a 
time t S ]0, oo[, where a: S M^ is position and p S M'^ denotes momentum. In the 
regime in which collisions among the particles can be neglected, the time evolution 
of the distribution function / is given by the Vlasov equation 

(1.1) dtf + v-V^f+{E + c-\xB)-Vpf^Q, V ^ 



where v is the relativistic velocity and c the speed of light. Here the mass and 
charge of the particles have been set to one. E — E{t,x) and B = B{t,x) denote 
the self-induced electric and magnetic fields, given by the Maxwell equations 

(1.2) V X B- c-^dtE = inc-^j, V -B = 0, 

(1.3) y X E + c^^dtB = 0, y-E :^ Anp. 
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The Vlasov and Maxwell equations are then coupled via the charge and current 
densities 

(1.4) p= fdp and j = / vfdp. 

Equations (ll.ip - (|1.4p are known as the relativistic Vlasov-Maxwell (RVM) sys- 
tem, which is essential in the study of dilute hot plasmas. Details and an abundant 
bibliography on this system can be found, for instance, in |6] . 

We further decompose the electric field into E = Ej^+Ex, where the longitudinal 
El and transverse Et components of the electric field satisfy, respectively 

(1.5) W X El^O and V • £'t = 0. 

If we now neglect the transverse component of the displacement current dtEr in the 
evolution equation (|1.2|) -the so-called Maxwell- Ampere equation-, then the RVM 
system reduces to 

(1.6) dtf + v-V^f+{EL+ET + c-^vxB)-Vpf^O, v^ ^ 



l + c-2|p|2 



coupled with 

(1.7) V X B - c~^dtEL = A-Kc-^j, V -B = 0, 

(1.8) V X Et + c^^dtB = 0, V -El = iirp, 

by means of (|1.4p . Equations (|1.4p - (ll.8l) are the RVD system. From the phys- 
ical point of view, the Darwin approximation is valid when the evolution of the 
electromagnetic field is 'slower' than the speed of light. 

In this paper we are concerned with the Cauchy problem for (|1.4p - (ll.8p . Global 
existence of weak solutions was shown in [5] for small initial data. The smallness 
assumption was later removed in |13j , where the existence and uniqueness of local 
in time classical solutions was also proved. In |15) , it is shown that solutions having 
the same regularity as the initial data (which is not the case in [13]), can be extended 
globally in time provided the initial data is small. At the present time, the existence 
of global in time classical solutions for arbitrary data remains unsolved. Here, we 
provide a constructive and somewhat simplified proof to the local in time existence 
and uniqueness result for classical solutions of the RVD system, and we show that 
the solutions can be extended for all times if the initial data are sufficiently small. 

The main difficulty when dealing with the RVD system has been to find an a- 
priori estimate on the transverse component of the electric field E^- In contrast 
to the RVM system, the component E-r does not contribute to the energy of the 
electromagnetic field, and thus the law for the conservation of the total energy does 
not provide any control on the L^-norm of Et- Indeed, the total energy of the RVD 
system reads 

c'Jl + c-2 \pff{t, x,p)dpdx + -^ / \\ELit,x)\^ + \B{t,x 



dx. 

Hence, by virtue of the underlying elliptic structure of the Darwin equations, duality 
type arguments and variational methods have previously been used to estimate Et- 
Here, we take advantage of the formulation of the RVD system in terms of the 
generalized variables -defined later on-, and we produce an L^-bound on p^/'^Et 
instead. This estimate is at the core of our results, and is given in Lemma [51 It is 
remarkable that by pursuing such an estimate we have obtained, 'almost for free'. 
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an L^-bound on dxEr as well. In contrast to the results cited above, the law for 
the conservation of the total energy is not used in our proofs at all. 

The structure of the paper is as follows. In Section [51 we present the scalar and 
vector potentials and introduce the generalized position and momentum variables. 
We then recast the Vlasov and Darwin equations in terms of the new variables, and 
treat them as uncoupled linear equations. A representation for the Darwin vector 
potential is given, and some standard a-priori bounds are obtained as well. Then, 
in Section [3] we couple both Vlasov and Darwin equations and introduce the RVD 
system in terms of the potentials. The estimates on the transverse component of 
the electric field and its space derivative are produced in Subsection 13.11 Finally, 
we study the Cauchy problem for the RVD system in Section |4l First, we produce 
the local in time existence result in Subsection 14.11 and then, in Subsection 14.21 
we extend local solutions globally in time under the smallness assumption on the 
Cauchy data. We conclude with an Appendix. 

We remark that the RVD is actually an hybrid system, since we are considering 
relativistic charged particles whose interaction with the electromagnetic field they 
induce is an order- (w/c)^ approximation |10|l9]. Yet, the RVD system is interesting 
in its own right, in particular for numerical simulations, since it contains an underly- 
ing elliptic feature while preserving a fully coupled magnetic field. This is in contrast 
to the more involved RVM system, whose hyperbolic structure yields both analyt- 
ical and numerical challenges. Also, the tools used here are likely to be adapted to 
the 'proper' physical system, which is (|1.4p - (|1.8|) with v = p (l — c~^p'^/2) instead. 

The following notations will be used in the paper. As usual, C^'"-{X; Y) denotes 
the space of functions / : AT ^ y of class C^ whose k-ih. derivatives are Holder 
continuous with exponent a G (0, 1). Cq{X\Y), resp. C^{X\Y), are the spaces of 
C'''(A; F)-functions with compact support, resp. bounded. W^'°"{X\Y) stands for 
the Sobolev space of L°°{X; y)-functions whose weak first order partial derivatives 
belong to L°°{X;Y). If I is an interval in R, then by g G C^ (/, C'^(A); F), we 
mean that g: IxX^Y,g = g{t, x), and for all i G /, g{t) G C''{X; Y) and the 
function 1 1— ^ g{t) G C^{X] Y) is of class C^ on /. For such a function, we sometimes 
write (by abuse of notations) g G C^{X\Y) to mean that git) G C^{X\Y) for all 
t € I. Similarly, the norm of g(t), say the L'^-norm ||5(i)||Lj, will sometimes be 
denoted by ||g||i9 . All other notations in the paper are standard, and the constants 
may change values from line to line. 

2. The Potential Representation 

From classical electrodynamics it is known that an electromagnetic field (-E, B) : 
M X M^ -> M^ X ]R3 that is a smooth solution of the MaxweU equations (fO |) - ([01) 
can be represented by a set of potentials ($, A) : (0, oo) x R'^ ^ M x M'^ according 
to the expressions 

(2.1) E{t,x) = -\7<i>{t,x)-c-^dtA{t,x), 

(2.2) B{t,x) = WxA{t,x). 

These relations can easily be obtained from the two homogeneous Maxwell equa- 
tions in (jl.2l) - (jl.3p . In particular (12. 2p follows from the vanishing divergence of 
the magnetic field, while (|2.ip follows after inserting (j2.2|) into the remaining ho- 
mogeneous equation. Since for any smooth scalar function A we have the identity 
V X VA = 0, it is clear that such potentials are not uniquely determined. We may 
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find another pair given by 

(2.3) ($',A') = ($-c-i9tA,A + VA) 

which also satisfies (J2.ip - (l2.2p . The two sets of potentials are fully equivalent, in 
the sense that they produce the same electric and magnetic fields. 

This lack of uniqueness allows to impose a condition on the potentials that 
ultimately determines their dynamical equations. Even after doing so, some arbi- 
trariness remains that can be avoided by imposing an additional restriction on A. 
The resulting restricted class is called a gauge, and all potentials within this class 
satisfy the same gauge condition. Commonly, the Lorentz gauge condition 

(2.4) V -A + cr^dt^^O, 
or the Coulomb gauge condition 

(2.5) V-A = 

is used. The former is relativistically covariant and leads to a class of scalar and 
vector potentials that satisfy wave equations. This is a natural choice when dealing 
with the RVM system. It was used in [3] to study the smoothing effect resulting 
from a coupling of a wave and transport equations. It was also used in [4] to produce 
an alternative proof of the celebrated result by Glassey and Strauss on the RVM 
system 0. On the other hand, the Coulomb gauge condition leads to scalar and 
vector potentials that satisfy a Poisson and a wave equation, respectively. As we 
shall see in Subsection l2.2l below. this is the correct choice to introduce the potential 
representation of the RVD system. In a way, both the Lorentz and Coulomb gauges 
can be seen as limit cases of a more general class known as the velocity gauge, in 
which the scalar potential propagates with an arbitrary speed [9|. 

2.1. The Vlasov Equation. We now introduce the generalized variables, which 
permit to rewrite the Vlasov equation (|l.ip in terms of the scalar and vector po- 
tentials in a very convenient way. The resulting transport equation is shown to 
be determined by an incompressible vector field irrespective of the gauge chosen. 
Thus, we can count on the usual a-priori estimates on the distribution function -see 
Lemma [2] below- no matter which gauge we decide to work in. 

To start with, let I C [0, oo[ such that G /. Assume that the pair {^,A) e 
C^ (/, C^ (R3); K X K^) is given, and so in view of (I2ll])-(l221) the electromagnetic field 
is given as well. Denote z := {x,p) and write {X, P){s) instead of (X, P){s, i, z) to 
ease notation. Then, by virtue of (|2.ip - (|2.2l) . the characteristic system associated 
to the Vlasov equation (|l.ip reads 

(2.6) X{s) - v{P{s)), 

(2.7) P{s) = [-W(^-c-^dsA + c-\x{W xA)]{s,X{s),P{s)). 

Hence, since 

A{s, X{s)) = [dsA + (« • V) A] {s, X{s)), 
the equation (j2.7p can be rewritten as 

(2.8) P{s) = \-c-^A - V$ + c-^v X (V X A) + c-^ {v ■ V) a] (s, X{s),P{s)). 

The structure of (|2.8p suggests that we can define a generalized momentum variable 
TT = p + c~^ A such that the above equation can be reduced to 

tl{s) = [-V$ + c-\ X (V X A) + c"i (v ■ V) A] (s, X{s),P{s)). 
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Here we have denoted n(s) = P{s) + c^^A{s,X{s)). On the other hand, the 
relativistic velocity written in terms of the generahzed momentum is 

/„ „N Tr-c~^A 

(2.9) VA = , 

Therefore, by using the elementary identity 

VAX{V xA) + {vA • V) A = vyVA\ 

we can reformulate the characteristic system (|2.6p - p.7p in terms of the generalized 
variables ^ = {x,tt) as 

(2.10) X(s,i,0 = VA{s,Xis,t,^),nis,t,0), 

(2.11) ri(s,i,0 = -[\/<P~c-'v\\/A']{s,X{s,t,OMs,t,0)- 

As usual, repeated index means summation. Now, standard results in the theory of 
first order ordinary differential equations imply that for every fixed < G / and ^ e R^ 
there exists a unique local solution S = {X,U){s,t,£_) of (|2.10p - (|2.1ip satisfying 
E{t,t,C) = ^; see [HI Chapters II and V]. Moreover, S e C^ (/ x / x M6;M'^). In 
turn, uniqueness implies that 

Z = (X, n - c"^v4)(s, t,x,TT- c-M) 

is the unique solution of (|2.6l) - (|2.7p with initial data Z{t,t,z) — (a;,7r — c^^A), so 
by having the characteristic curves in the generalized phase space we can recover 
the characteristic curves in the usual phase space. 

As the following lemma shows, the field resulting in the right-hand side of the 
system of equations (|2.10p - ()2.1ip is an incompressible vector field: 

Lemma 1. For va given by h2. ff|) . we have 

V^ • WA + V^ • (-V$ + c-^v\VA') = 0. 

Proof. Since trivially V^ • V$ = 0, the result is a consequence of the elementary 
relation 

D 

As a result, solutions of the characteristic system (|2.10p - (|2.1ip satisfy the volume 
preserving property. Specifically, for any fixed s, t G /, the map E!(s, i, •) : R^ — ;> R^ 
is a C^-diffeomorphism with inverse S~^(s,i,^) — 'B(t,s,^) and Jacobian determi- 
nant; see [8j Corollary V.3.1] 

detJ..,.(0 = ^^ = l. 

These properties of the characteristic flow lead to the following result: 

Lemma 2. Let ($,A) e C(/, C^(R'^);R x R^*) be given in some gauge and let va 
be given by h2.9\) . Assume that V<I> and V^', i = 1, 2, 3 are bounded on J xM.^ for 
every compact subinterval J <Z I . Let /o € C^ (R^; R) and denote by'E. = (AT, 11) the 
characteristic flow solving I12.10\) - ^2.11\) . Then, the function f{t,^) — /o(S(0,t,^)) 
defined on / x R^ is the unique C^ solution of the Cauchy problem for 

(2.12) dtf + VA ■ yj - [V$ - c-'v\VA'] ■ V,/ = 0. 
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Moreover, if fo>0 then / > 0. Also, for t <E I we have that 

suppf{t) = S(i,0, suppfo), 
and for each 1 < q < oo, t ^ I we have 

II/WIIlj,,-II/oIL.„. 

Conversely, if f is a C^ solution of the Cauchy problem for \2.12]) . then f is 
constant along each solution of the characteristic system i2.10\) - ![2.11\) . 

Remark 1. In addition, if ($, A) (i) e C2^"(]R3;R x R^), < a < I, t E I, and 
/o e C^'" (R6;K)^ then the unique C^ solution f{t,^) = /o(S(0,i,O) of the Cauchy 
problem for ([^1^ satisfies f{t) e C1'"(R6;M) for every t e /. 

Proof of Lemnia[2j In view of Lemmall] the proof follows by the standard Cauchy's 
method of characteristics; see [SI Chapter VI]. In particular, the properties of / are 
a direct consequence of the properties of the characteristic flow discussed above. D 
We point out that (|2.12p is the proper Hamiltonian representation of the Vlasov 
equation (|l.ip in terms of the potentials, since the characteristic equations (|2.10p - 
(|2.11l) are Hamilton's equations for the Hamiltonian 

(2.13) n{t, X, tt) = c^y^l + c-2 \n - c-^A{t, x)\^ + $(i, x) 

of a relativistic charged particle under the influence of an electromagnetic field of 
potentials ($, A). As before, in (|2.13p the charge and mass of the particle have 
been set to one. 

2.2. The Darwin Potentials. To determine the dynamical equations satisfied 
by the potentials we shall impose the Coulomb gauge condition, since it leads to 
the Darwin approximation of the Maxwell equations and ultimately to the RVD 
system. Throughout this section, unless we specify otherwise, we assume that both 
the charge and current densities p and j are smooth and given, and they satisfy the 
continuity equation 

(2.14) dtp + V-j = 0. 

Formally, if we substitute the electric and magnetic fields in (|2.ip - (l2.2p into the 
non- homogeneous Maxwell equations in (|1.2p - (|1.3p . we find that $ and A satisfy 

(2.15) A$ = -inp-c-^dtiV -A), 

(2.16) AA - c-^d^A = -c-Hnj + V (V ■A + c-^dt^) . 
Therefore, in the Coulomb gauge (|2.5|) . the potentials satisfy 

(2.17) A$ = -47rp, 

(2.18) AA-c-^d^A = -c-Hnj + c-^Vdt^. 

On the other hand, any smooth solution ($, A) of the above system that satisfies the 
Coulomb gauge condition initially, will continue to do so for all times, and therefore 
the induced electromagnetic field will solve (|1.2p - (ll.3p . Indeed, if ($, A) is a smooth 
solution of (PT7l) - (P?T51) that satisfies V • A\^^^^ = and dt{V ■ A)\^^^ = 0, then 
gc = V • A is the solution of 

Age - c'^d^gc = -4^c"i (V • j + dtp) = 0, 
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and the claim follows. Hence, the system of equations (|2.17p - p.l8p complemented 
with (|2.5|) is fully equivalent to the set of Maxwell equations ()1.2|) - (|1.3|) . 

We define the Darwin approximation of the Maxwell equations as the quasi-static 
limit of the system ^JT^-^A^: 

Definition 1. Let {p,j) : / x M'^ — > M x M'^ be given and satisfy the continuity 
equation \2.14^ . The set of potentials (<&, A) is called a classical solution of the 
Darwin equations if ^ e C^/, C2(M3); M), A e C{I,C^{M.^);R^) and, on I x M.^ , 

(2.19) A$ = -47rp, 

(2.20) AA = -c-^4:TTJ + c-^Vdt<^. 

The system (|2.19p - (|2.20p has the following explicit solution, as proved below: 

Definition 2. For the charge and current densities {p,j) : / x R'^ -^ M x R'^ we 
formally define the set of Darwin potentials ($/), A^)) : / x M"^ ^^ R"^ x R'^ by 

(2.21) <^D{t,x) = / p{t,y) '^y 



\y-xV 

dy 



AD{t,x) ^ — [id + uj®u]j{t,y) 
2c Jk3 



d3 



(,2.22j ^^^y^,^, . |^.„ , „^„jj^„,y^ . 

2c Jk3 \y - x\ 

where lo — [y — x) / \y — x\ and id denotes the identity matrix. 

Lemma 3. Let p e Ci(/, C^(R3);R) and j £ C(/, Co'"(R3);R3), < a < 1, 6e 
given -they do not need to satisfy the continuity equation {2.1^^ -. Define the field 

(2.23) ¥j{t,x)=j{t,x) + ^V f V-Jit,y)-^^, tel, xe 

47r Ja3 \y - A 

Then the following holds: 

(a): The scalar potential $c is the unique solution in C"'^(/, C^'"(M^); R) of 

(2.24) A$(i,a;) = -47rp(i,x), lim $(t,a;)=0. 

|a;|— >-oo 

It satisfies 

V^D{t,x)^ f p(i,2/)^^. 

JR3 \y — x\ 

(b): Fj e C(/,C1'"(R3);R3). It satisfies V -fj = (i.e., Pj is the transverse 

component of the current density j), and Vj{x) = 0{\x\ ) for \x\ — >■ oo. 
(c): The vector potential Aq is the unique solution in C(/, C'^'"(R^); R^) of 

(2.25) Ayl(t,a;) = -47rc^¥j(i,a;), lim |A(t,a;)|=0. 

|2;|— >oo 

It satisfies 

(2.26) dxAD{t,x)^—\ {uj®i~i®uj + [iuj®uj-id]{j-uj)}-^ ^, 



2c J ^3 \y — x\'^ 

with j = j(t, y). In particular, 

V-A£)(i,x)=0 and \7 y. Aoit^x) ^ - j uj x j{t,y) ^-^. 

c Jr3 \y — x\ 
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Corollary 1. If p and j, as given in Lemma \^ satisfy the continuity equation 
[2l4\ l, then 

Pj{t,x)^j{t,x)-^Vdt<^D{t,x), tel, xeM.^ 
An 

and thus the Darwin potentials li2.21]) - ^2.2^) are the unique classical solution of the 

Darwin equations \2.19j) - P2.20\) . 

Proof of Lemma 131 Without loss of generality we omit the time dependence. 

The proof of (a) is a standard result for the Poisson equation. Existence (in 
a much weaker sense) can be found, for instance, in 11, Theorem 6.21] while the 
regularity of the solution is given in [Hi Theorem 10.3]. Uniqueness is known as 
Liouville's theorem \\T, Theorem 7 Section 4.2]. 

To prove (b), notice that V • j e C^(M'^; M). Hence, as in (a), the integral in the 
right-hand side of (I2.23P is the C^'"-solution of the Poisson equation Au — —AttW-j, 
limu|_j.oQ m(x) = 0. That Vu(x) = 0(|a;|~ ) for \x\ — ^ cxd is well known, which in 
turn provides the decay for Pj, since j has compact support. Moreover, 

V • ¥j{x) = V • j{x) + ^A f V • j(y)^^ = 0. 
47r V \y - x\ 

As for (c), we first prove the following lemma: 
Lemma 4. The Darwin potential Ad in \2. 22]) has the equivalent representation 

(2.27) Ad{x) ^-f j{y)-^ + ^ [ ^- Mr—^.dy. 

c Jk^ \v-x\ 2cJsz \v-x\ 

Proof. The current density j has compact support, so standard arguments can 
show that the right-hand side (RHS) of the above expression is well defined. The 
divergence theorem then yields, 

WS = - f jiy)-^ -^ f ijiy) ■ V)udy 
c Jr3 \y - x\ 2c 7r3 

iiv)- 7;[j{v)-^{i{v)-^)] 



c Jk^ I 2 ] \y-x\ 

^ ^ [3{y)+Lo{j{y)-u)] '^y 



2c 7r3 \y - x\ 

which is precisely the Darwin potential Ad in (j2.22|) . The use of the divergence 
theorem is justified by the following standard argument: remove a small ball about 
a; S M^ in the domain of integration so we can avoid the singularity at y — x, then 
use the divergence theorem and note that the boundary term corresponding to the 
small ball vanishes as its radii tends to 0. D 

We shall now deduce by direct computation from (I2.27p . the Poisson equation 
given by (I2.25p . To this end, we first recall that just as in part (a), 

(2.28) A |i/,(.)^| = --.(.). 

I c Jr3 \y-x\} c 

The integral in curly brackets is in C'^'"(R^;R'^), due to the regularity of j. Next, 
we show that the following equality holds in the sense of distribution, 

dy 



(2.29) d,, / V • Jiy)^"dy = - / V • j(y) [6,k - cu'l^ 



I?/- 
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Let r = \y — x\ > 0. First, note that ^k^w' = —r^^ \5ik — cj'o;'^] , and that the 
integral on the right-hand side of (|2.29p is wcU defined for almost all x € R"^, since 
V-j e C^(R3, R) and the kernel is bounded from above by r^i. For (j) e C^(R3; R), 
the function {x,y) i~> dx^(t>{x)V ■ j{y)i^^ is integrable on R'^ x R^. Hence, we can 
use Fubini's theorem to find that 

d,,^{x)\ f W-j{yWdy]dx = f \ f {dxM^))w'dx]w-j{y)dy 

IJR3 J JR3 LJR3 J 

(t){x)dx^uj''dx\V ■ j{y)dy, 



where the second equality is justified by a standard limiting process and integrations 
by parts, similar to the argument at the end of the proof of LemmalU Then, another 
use of Fubini's theorem yields (|2.29l) in the sense of distribution, as claimed. 

Actually, the equality in (|2.29p holds in the classical sense. By the standard 
theory of the Poisson equation, the right-hand side of (j2.29p is a function in 
C^'"(R^;R); see [HI Theorem 10.3]. Therefore, in view of the theorem for the 
equivalence of classical and distributional derivatives, the integral in curly brackets 
on the left-hand side of (|2:29)) is in C3'"(R3;R); see [HI Theorem 6.10]. 

Now, since dx^T~^ — r~^uj'' and w*^ [5ik — w^w'^] = 0, we have 



Therefore, similar arguments to those used above yield 



^, ,fci dy 



I Jks \y 

(2.30) = -2/ V.j(y)^^=-2a.. / V.j(y)^^. 

Jk.^ \y-x\ 7r3 \v~x\ 

Hence, since A = V • V, we can combine (|2.29p and (|2.30l) to find that 

(2.31) a(1 / V.j(y)f^dyl = -iy / V • ,(2/)/^ 



ZcJrs \y~x\ ] c Jts.3 \y-x\ 

Then, we add (l^?^ and ([^31]) to conclude that AA^ = -Anc-^Fj holds on R^, 
and so Ajj is a C^'°' solution of (12.251) . This solution is unique in view of the 
Liouville's theorem [12, Theorem 7 Section 4.2]. 

The representation (j2.26p of dxA can be proved as follows. Since ja is regular 
enough, we shift the x-variable into the argument of ja and differentiate (|2.22p 
under the integral. Then, we move the derivative to the kernel of (I2.22p helped by 
the same standard argument at the end of the proof of Lemma 21 In doing so, we 
notice that for r > 0, the imk-th entry of dxJC is 

which leads to (|2.26p . Finally, it is not difficult to check that 



1 



\/ ■Ad= Tra.ce{dxAD) = and (V x Ad)' - - [dxAo - {OxAd 



\T\ 



2 



'kl ■ 



where {dxAo)^ denotes the transpose of dxAo, and i,k,l e {1,2,3} are given 
according to the cyclic index-permutation. D 
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It is easy to check that the Darwin equations in Definition [T] are formally equiv- 
alent to the equations ()1.5p and ()1.7|) - (|1.8p given in the Introduction. To see this 
let us define 

We have to show that {El,Et,B) formally solves (jl.7p - (|1.8l) provided that the 
charge and current densities satisfy the continuity equation. Clearly V ■ B = 0, and 
since V • Ajj = 0, we have 

V X B = V (V • Ad) - AAd = inc-'^j - c-^a^V^D = 47rc-^j + cr^dtEL, 

which is (|1.7p . Easy computations yield (|1.5p and (|1.8p . and the claim follows. 

We conclude this section with a-priori estimates on the potentials and their space 
derivatives. For simplicity and without loss of generality, we shall neglect the time 
dependence. 

Lemma 5. For 1 < m < Z set vq = 3/(3 — m) and let r < vq < s. Then there 
exists a positive constant C = C{m, r, s) such that for any ^ G L*" n iy'*(M^; R) 



nvh '' 



\y 



< C(m,r,s) 11*11^7^ 11*11^, , where X-^ ''^''^ 



1 — r/s 

In particular, C{m, 1, oo) = 3 (47r/m)™ / (3 — m). 

Proof, cf. 13, Lemma 2.7]. D 



Lemma 6. For p and j as given in Lemma\^ the Darwin vector potential J 
satisfy the estimates: 

(2.32) \\Aoh^<C\\j\f/l\\j\\'/l and \\d^Aoh^ < C \\j\\'l^ \\jf^^ 

Moreover, for any < h < R we have 



diAnL^<C 



iDlli^ -o R-^\\j\\^^+h\\d.,j\\^^ + il + \niR/h))\\j\\^^ 

where C > is independent of h, R, p and j. In particular, 

(2.33) \\dlAn\\^^ < C [||j||^, + (l + \\j\\^^) (l + ln+ \\d,jh^) 

The same estimates hold for the scalar potential $c with j replaced by p. 

Proof. The estimates corresponding to $£i are well known from the study of the 
Vlasov-Poisson system. These results can be found, for instance, in [14j Lemma PI] 
and [l, Propositions 1 and 2]. Here, we shall produce the estimates for the vector 
potential A]j only. The proof is actually rather similar. 

Let IC{y,x) — \y — x\~ [id + a; ig) w]. Clearly, |/C(y,a;)| < C |?/ — a;|~ . Then, the 
estimates in (|2.32p are a straightforward consequence of Lemma [5l To produce the 
estimates for the second derivatives, consider 

\y-x\ 
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Here we have introduced the notation dk = dx^, fc = 1, 2, 3; see Lemma[3l^c) for the 
matrix representation of the integrand of dxA. Now, the integral /i can in turn be 
spht into three integrals, each one essentially the same as the integral corresponding 
to didk^D- Thus, /i satisfies the expected estimates, as didk^o does. Therefore, 
we are led to estimate l2- To this end, we set r — \y — x\, and for r > we denote 
hy T^iiv-x) 



a, 



w^w^w*^ 



= — [(^miW^Cj'' + SuUj''uj"^ + Sklio'-Uj"^ - Scj'w'^ILj'cj 



This kernel is too singular to use Lemma [SJ However, since y^y^y"'' \y\ is homo- 
geneous of degree —2, for every < i?i < i?2 we have 



/ 



^Tki{y)dv 



fi> I i5 "''^y 



lRi<\y\<R2 J\y\=R2 -"-2 \y 

Thus, for any /i > 0, we can rewrite I2 as 



yl yiykym ^^ 

\y\=Ri ^1 \yt 



0. 



h = / T2i{y-x)j'^{y)dy + j'^(x) uj' lo^ J uj"" duj 

>h J \uj\ — l 

TTu{y-x)[r{y)-r{x)]dy. 

\y~x\<h 

The singularity in the last integral at r = is now avoided by the difference 
j"^{y) — j"^{x). Indeed, for < /i < i? we produce 



h < C' 



IIjIIl= 



\dxj\ 



dy 



h<\y~x\<R \y- X\ 

dy 



y-x\>R \y-x\ 



\y-x\<h \y-x\" 



W)\ 



< c 



HR/h)\\j\ 



R-'\\j\\L^+h\\dxjh^ + \\jh. 



This yields the first estimate on ||9^^d| 



Then, by setting R = 1 and letting 



l|c^a;j|lL=° ^^ W'^xJWl^ > 1, Otherwise h = 1, the estimate (|2.33p follows as well. 



This completes the proof of the lemma. 

3. The RVD System 



D 



If we now combine (|2.12p and (|2.2ip - (|2.22l) by means of (|1.4p . then we obtain 
the following equivalent representation of the RVD system: 

(3.1) dtf + VA ■ Vxf - [V$ - c-\\VA'] ■ Vp/ = 0, 

coupled with 

dpdy 



(3.2) 

(3.3) 
where 

(3.4) 



$(i,x) 
A{t,x) 



fit,y,p) 



\y-x\' 



1 



VA 



[ld + uj ® uj]vAf{t,y,p) 



p- c^^A 



dpdy 

\y-x\' 



^l + c-^\p-c-^A\' 
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For the sake of notation, we have written p instead of n when referring to the 
generahzed momentum variable. We wiU continue to do so for the rest of the paper. 
We shaU also set c = 1 for the speed of light. The goal is to prove that a small 
enough Cauchy datum launches a unique classical solution of the system (|3.1[) - (I3.4|) 
globally in time. We shall prove this in Subsections 14.11 and 14.21 below, but first 
we center our attention on p.3|) . If / is given, then p.3|) is a nonlinear integral 
equation of unknown A. 

Lemma 7. Fix t e I and let f{t) e Cg'"(]R^;M), < a < 1, be given. Then, there 
exists an A{t) € Ct n C2'"(M3;M3) satisfying (E3)-(3^- 

Proof. Without loss of generality we shall omit the dependence in time. Let C be 
a constant that may depend on /, to be fixed later on. Define the set 



V, 



{AeCfc(K3;K-'^) : \\A\\^^ < c} 



c 



First, we show that there exists an Aqo G T^q which solves (|3.3p - p.4p . To this end. 



denote the kernel ]C{x,y) = \y — x\ 
mapping ^ i~> T[A] defined by 



T[A]{x) = 



1 



K^ix,y)vAf{y,p)dpdy, va 



uj] and let A G Vq. Consider the 
p- A 



l + \p-AY 



We claim that T[A\ e V^. Indeed, let {IC)ij{x,y) be the zj-entry of IC{x,y). For 
some Ml, U2 and U3 on the line segment between x and z, the mean value theorem 
implies 



i^^)^Ax,y)- {K.)iAz,y) 



\y-ui 

Hence, since \va\ < 1, a use of Lemma 12] produces 

1 




y~u3\ 



\T[A]ix)-T[A]{z)\ 



(3.5) 



< 



< 



< 



2 Jr3 
C\x 

C{f)\x^z\ 



\IC{x,y)~IC{z,y)\p{y)dy 
dy 



p{y): 



\y 



Thus, T[A\ is a continuous vector valued function. Also, by Lemma [6] 



(3.6) 



|T[A]||^^<3/2(7r/2)V-'||p||^;;^||p||i(^^C 



2/3 



1 1/3 _ 



Therefore, T\A\ G Vq as claimed. 

We now show that T has a fixed point A^a G T^c- ^Y virtue of the Schauder 
fixed point theorem [12j Theorem 3 Section 9.1], it suffices to show that T is a 
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continuous mapping and that the closure of the image of T is compact in "Dq. To 
show the continuity of T, we see that if Ak — > A in !?(=;, then by Lemma [6] 

\T[A^]{x)-T[A]{x)\ < cf / bA.-«A|/(y,p)i*'^^ 



< C{f)\\Ak~A\\^^. 



To show that TVq C Vq is compact, we first notice that for A £1), 



(3.7) \T[A]{x)\<\\p\\^^ I ^^<C{f)- ^ 



supp/ \x-y\ 1 + \x\ 

Now consider the sequence {i?„} C TV^. Let i? > be fixed. By (j3.6p and (|3.5p . 
the restriction 

{^n}\[xeS.^:\x\<R} 

is clearly bounded and equicontinuous. Then, by Arzela-Ascoli and a standard 
diagonal argument we can find a subsequence {Bn^} and a continuous, bounded 
limit vector field B such that {_B„^. } — >■ B uniformly on compact sets, and in 
particular pointwise. Clearly, ||i3||^oo < C, and since {-Bn^} satisfies the estimate 
p.7p . so does B. We only need to show that the convergence {Bn^ } — > S is uniform. 
To this end, let e > 0. Choose R > such that the right-hand side of (13. 7p is less 
than e/2 for |.-e| > R. Then, for all k we have |i?„j.(a;) — B{x)\ < e for |a;| > R, and 
we can find a, kg = fco(-R, e) such that for all k > ko 

sup |B„fe(a;) - B{x)\ < e. 

\x\<R 

This proves uniform convergence. Hence, all the hypotheses for the Schauder fixed 
point theorem are fulfilled, and thus T has a fixed point Aoo in V^. 

Next, we have to show that Aoo has the required regularity. To this end, define 
VA^ and then ja^ according to (j3.4|) and (II. 4|) . respectively. The vector field Aoo 
has the form of a Darwin potential (|2.22p with current density ja^ € Co(K."^;M''). 
Clearly, the kernel of (|2.22p satisfies \IC{x,y)\ < C \y — x\~ and the derivative 
estimate \dxJC{x,y)\ < C \y — x\ . Hence, we can use the standard theory for the 
Poisson equation to find that Aoo G C^(K^; K"^); see, for instance, [U Lemma 4.1] or 
[TTl Theorem 10.2 (in)]. But such a regularity of ^oo imphes that ja^ € Cg (R^; M^). 
Thus, we also have ja^ € C^(]R3;]R3)^ < a < 1 and so ^oo e C'2'"(M3;M3)^ ^^g 
desired. For the latter implication see, for instance, jll) Theorem 10.3]. D 

Remark 2. If wc consider the time dependence in Lemma [71 and assume that / is 
C^ with respect to t e /, then A is also C^ in i G / as a consequence of the Implicit 
Function Theorem in Banach spaces; see |16] . 

3.1. Estimates on dtA and its space derivative. We now turn to the estimates 
on the time derivative of the vector potential (i.e., the transverse component of the 
electric field), and its space derivatives. 

Throughout this section, we shall assume that the triplet (/, $, A) satisfies (|3.ip - 
(|3.4p on / X R-^ X M'^ with f{t) having compact support on M'^ x M'^. For / as given, 
define t M- Z(i) by 

(3.8) Zit)=sup{\ix,p)\ ■.30<s<t:fis,x,p)^0}. 
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The function Z{t) is a non-decreasing function of t, which by the compact support 
of / is bounded on any finite subinterval of /. The following lemma is essential to 
our results: 

Lemmas. Let f e C^/, Co'"(K^);K) and{^,A) € C^{I,C^^°'{M.^)-R^ xR^), with 
/ > and < a < 1, satisfy i3.1] )- (K4^ . Define p and Z{t) according to ^1-4^ and 
\3.8]) . respectively. There exists a positive C(t) = C{Z{t), \\f{t)\\j^o<D ) such that 

\\dtd^A{t)\\^, + p'/^{t)dtA{t) <C{t), tel. 

Remark 3. For t G I, Lemma [S] and the assumption on the support of / imply 
\p- A\ < C{Z{t), ||/(i)|lroo ) < oo and so \va\ < 1 strictly on supp/(i). 

Proof of Lemma [51 For va as given in (|3.4p define the current density 

JAit,x)^ / VAf{t,X,p)dp. 



By Lemma [DJc), the components A'', i — 1,2,3 oi the vector potential satisfy 
(3.9) AA'it,x)^-A7TJ\it,x)-d,J ^.j^lt,y)-^. 

Jr^ \y-x\ 

Take the partial time derivative on both sides of the above equation and multiply 
by dtA^. After integration by parts, dropping the An and using the definition oi ja, 

IVdtA'l^ {t,x)dx ^ f f dtA\t,x)dt{v\f){t,x,p)dpdx 

dxdy 



(3.10) - / / dtd,^A\t,x)V ■dtJAit.y) 



\y-x\ 



Note that the boundary terms vanish. Indeed, since / has a compact support, 
so does JA and the boundary term corresponding to the first term on the right- 
hand side of the above equation is zero. On the other hand, it follows by standard 
arguments that dtA^{x) has at least a decay Odxp ) and VdtA^{x) = 0(|a;|~ ). 
Moreover, the integral I{x) on the right-hand side of (|3.9p has a decay 0(|a;|~ ) and 
so does dtl{x). Therefore, dtA^^dtA'{x) = 0{\x\~^) and dtA'dJix) = 0{\x\^^), 
which suffice for the boundary terms to vanish. 

Now we add the equations p.lOp for each component of A. We find 

\d^dtA\' = / f f{dtA-dtVA)+ f f {dtA-VA)dtf 



-dt [ / - (V • A) (V • duA) 

JR3 JR3 r 



(3.11) = h+h + I: 



3- 



But /3 = since the vector potential satisfies the Coulomb gauge condition; see 
Lemma [3Uc). Also, by using the representation of the derivatives of the velocity 
given in the Appendix, the integral /i can be written as 

h = - f [ -^=L=(\dtA\'-\vA-dtAf 
jm3 Jm^ yi + 5 ^ 
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where we have denoted g — \p — A\. We shaU also denote K = — V<i> + v\WA\ 
Hence, after sending /i to the left-hand side of p.ll|) . and by using the Vlasov 
equation (|3.ip in I2, we find that 

/ 



\d,dtA\' 



" V^l + 9' 



{\dtA\ 



' \vA ■ dtAf 



(3.12) 



{VA ■ dtA) [V, • (vAf) + Vp • (Kf)] 
f f fdtA' {vA ■ W^v\) + f f fv\ {vA ■ W^dtA') 
+ / / fdtA' {K ■ S/pv\) . 



Notice the integration by parts and the use of the product rule in the last equality. 
We claim that the left-hand side of the above equality has a positive lower bound 
for every time t. Indeed, we have that 



\dtAf - \vA ■ dtA\^ > \dtAf - \va\^ \dtAf 



1 



\vAn\dtA\\ 



Also, by Remark [3] there exists a .gmax(^) = 5'max(^(^), 11/(^)11 roc ) < oo so that 

(3.13) i^tf . — ^ > '-^ ^ a^n > 0. 



Therefore, the left-hand side of (I3.12p satisfies 



(3.14) 



LHS > C^Ut) i \\dxdtA{t)\\l, + p"\t)dtA{t) 



On the other hand, the known bounds on the derivatives of the potentials given in 
Lemma m lead to 

\\d.VA{t)\\L^ +\\dpVAmL- +\\KmL^ <C{t)^C{Z{t),\\f{t)\\^^)- 

see the Appendix for an explicit representation of the derivatives of the velocity. 
Hence, after a use of the Cauchy-Schwarz inequality and again the use of the com- 
pact support of /, the right-hand side of p. 121) can be estimated as 



(3.15) RHS < C{t) i\\d,dtA{t)\\^, + p^/^{t)dtA{t) 

Finally, since (a + hf <2 {a^ + P), the resuh follows from (|XTi)) - (|XT5)) . 
Lemma 9. Under the assumptions of Leninia\^ we have that 



D 



\dtA{t)\\ 



< C 



\pm]!'Mt)f4{i + \\pm%'\\pm'4 



\\pm'i!'\\pmi^' p'/'imA{t) 



Ll 



\L\. Wl^yiWL 

Corollary 2. Under the assumptions oj Lemma\^ we have that 

\\dtA{t)\\^^<C{t), tel, 

for some positive C{t) — C{Z(t), ||/(i)||^oo )■ 



tei. 
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Proof of Lemma [9l Consider the integral representation p.Sp for the vector poten- 
tial, and take the partial time derivative on both sides of this equation. Denoting 
the kernel by /C(a;, y) and dropping the multiple 1/2, we have 



dtA{t,x) = IC{x,y)[vAdtf + dtVAf]{t,y,p)dpdy 

= h+h- 

Set K = — V<i> + w^VA*. A use of the Vlasov equation yields 

h ^ - I f JC{x,y)vAVy{vAf)- f f JCix,y)vAVp-{Kf) 



fvA-[dylC{x,y)vA+IC{x,y)dyVA]+ / fK-K{x,y)dpVA- 

R3 JR3 JR3 

Therefore, since \va\ < 1, also \dxJC{x^y)\ < C \y — x\ (see Lemma [3jc)), and 
|9a;WA| < C l^a;^! and |9pWA| < C (see Appendix), we have 

h < cU p{t,y)-^ + (\\dMt)\\,^ + \\dxAmL^) [ pit,y)T;^] 

[jR3 \y-x\ ^ " " ^ jR3 \y^x\) 

(3.i6)< ciip(i)iiif iip(i)iii§(i + iip(i)iii('iip(i)iii(?), 

where in the last inequality we used the estimates from Lemmas [5] and |6l 

On the other hand, since |9tWA| < C* |9tv4| (see Appendix), the integral I2 can 
be estimated as 

h < C j p{t,y)\dtA{t,y)\. '^y 



\y~x\ 

< c\\p{t)d,Am]!,^ \\p{mA{t)\\%^ , 

where the second inequality is a consequence of Lemma [5] Hence, the Cauchy- 
Schwarz inequality and a direct estimate lead to 

(3.17) h < c\\pmi{'\\pit)\\'/i p'^'imAit) 

The lemma then follows from p.l6p and (|3.17l) . D 

Lemma 10. Under the assumptions of Leninia\^ we also have 

\\dtdxA{t)\\^^ < C{t) {\\dtf{t)\\^^^ + \\p{t)\\l{^ Mt)\f4 \\9tA{t)h^) ,tel. 

for some C{t) ^ C{Z{t)). 

Corollary 3. Under the assumptions of Lemma\^ we have that 

\\dtdxA{t)\\^^<C{t), tel, 
for some positive C{t) = C{Z{t), \\f{t)\\^^ , Wdtfm^^ ). 
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Proof of Lemma [TOl Consider 

dtdxA{t,x) = / / dxK:{x,y)[vAdtf + dtVAf]{t,y,p)dpdy 

(3.18) - h+h, 

where \dxJC(x, y)\ < C \y — x\ . Hence, by Lemma [5l we obtain 

h < C{Z{t)) \\dtI{t)\\L^ and h<C \\dtA{t)\\^^ IIpWH^^' \W)\fI^ , 
where |9tWyi| < C \dtA\ has been used. The result readily follows. D 

4. The Cauchy Problem for the RVD System 

A noticeable advantage of writing the RVD system in terms of the generalized 
variables and potentials is that it resembles, to some extent, the well-known Vlasov- 
Poisson (VP) system. Actually, the latter can be formally obtained from (|3.1[) - (l3.4p 
by letting c — )■ oo, so that terms involving the vector potential are no longer present. 
This resemblance allows to adapt previous techniques used for the VP system to 
the Darwin case. Below, the proofs we present are in the same vein as those given 
in [13] for the VP system. Obviously, several non-trivial difficulties arise due to the 
inclusion of the vector potential in the system equations, not present in the Poisson 
case. Incidentally, we expect that a global in time existence result to the relativistic 
Vlasov-Poisson system for unrestricted Cauchy data, which is still unsolved, will 
lead to an analogous result for the RVD system. 

4.1. Local Solutions. In this section we shall produce a local in time existence 
and uniqueness result for classical solutions of the RVD system. 

Definition 3. Let /o be given. We call f a classical solution of the RVD system 
if f e C^{I X R6;R); it induces the potentials ($,A) e C^{I,C'^{R^)]R x R^) ^^ 
i ti^.i^p -i TO)) .' for every compact interval J <Z I the fields V<1> and v\V A^ are hounded 
on J xR? and J x R'^ x R^ respectively; and the triplet (/, $, A) satisfies the system 
\S.l\) - [3^^ on / X R'^ X R'^. Moreover, we say that f is a classical solution of the 
Cauchy problem if /|j^q = fo- 

Theorem 1. Let /o G Co'"(K^;R), < a < 1, /o > 0. For some T > 0, there 
exists a unique classical solution f on [0,T[ of the Cauchy problem for the RVD 
system. Moreover, for each < t < T, the function fit) is in C"'^'"(K^; R), it is 
non-negative and has compact support. In addition, if T > is the life span of f , 
then 

sup{|p| -.30 < t <T,xeM.^ : f{t,x,p) 7^ 0} < oo 
implies that the solution is global in time, i.e., T — 00. 

Uniqueness . Consider two solutions (/i, $1, Ai) and (/2, $2, ^2) of the RVD sys- 
tem as given by Theorem [TJ The Vlasov equation yields 

dt{fl-f2f + VA,-^Ah~f2f+Kl-Wp{fi-f2f 

- 2 (/i - /2) [{VA, -VA,)- V,/2 + (i^2 - i^l) • Vp/2] 

where Ki — — V$i -|- v\ V A\ and analogously for K2. In view of the compact 
support of the solutions /i and /2, we set i? > such that 

supp/i(i) U supp/2(t) cBrxBr, te [0, T] C [0, T[. 
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Let Q{t) = \\fi{t) - f2it)\\l2 ■ Since /2 G CH[0,T[xR6;K) has compact support, 
||Va;/2(t)||ioc +\\'^pf2{t)\\rJ <Cr. Also, we have \va, - vaJ < CjAi -^2! and, 
by Lemmaini ||9a;(Ai,742)||£^oo < Cr. Then, it is not difRcult to check that 

(4.1) + ||Ai(i) - A2(i)||i.(5^) + iia.Ai(i) - a.A2(t)||i.(s^) 

From the Poisson equation satisfied by the scalar potentials we deduce 

/ \d^^,{t,x)\'dx^Ai, [ [ p,{t,x)p^{t,y)^^^, 

Je3 Js.3 Jr3 \y - x\ 

and analogously for $2. Linearity, the Hardy-Littlewood-Sobolev inequality and 
Jensen's inequality yield 

(4.2) |!9.$i(t) - 9.*2Wlli2 < C \\pi{t) - P2WL./5 < CflQ'/'W. 

On the other hand, in order to estimate the terms involving the vector potential, 
we proceed as follows. Define fx = A/i + (1 — A) /2 for < A < 1. Clearly /a > 
has compact support and satisfies d\f\ = /i — /2- Let A\ be the Darwin vector 
potential induced by fx ■ In view of Lemma [3] we have 

AAxit,x)^~A7: [ VAjxit,x,p)dp^V [ [ V . (vaJx) {t,y,p)^^^, 

J Br J Ba J Br \y ~ ^\ 

and V ■ Ax — Q- Notice that Ai (resp. A2) solves the above equation when A = 1 
(resp. A = 0). By virtue of Remark [2] (where t is replaced by A), we can use the 
arguments in the proof of Lemma [5] to find 

/ \dxd^Ax\^ + [ [ J": (\dxAxf - \vA, ■ dxAxf 

JK3 J Br J Br V^ +.9A \ 

(4.3) ^ f f i9^Ax-VA,)dxfx. 

J Br J Br 

The analogous expression in Lemma [5] is the first equality in (I3.12p with the right- 
hand side replaced by the expression of I2 in (|3.11l) . Note the integration over 
Br X Br in view of the compact support of dxfx- Hence, since \vax\ < 1 by 
Remark [3l we can use again the arguments in Lemma [8] and the Cauchy-Schwarz 
inequality on the right-hand side of (|4.3p to obtain 

Wdxd^Aximl, + pl/\t)dxAx{t) ]^<CR\\dxAx{t)h,^s^^\\dxfxmLi^, 
which implies that 

(4.4) \\dxd.Ax{t)\\l, < Cr ||aAAA(i)IL.(B„) Q'^'it) 

for some Cr > and all < A < 1. Poincare's inequality and (|4.4I) then yield 

WdxAximliiBR) < CRWdAAximl. < CR\\dxAxm^.^^sR)Q'^^it), 
and thus, 

(4.5) \\dxAxm^,js^^<CRQ'/^{t) 

for all < A < 1. Inserting (H3)) into dO]) . we also have for all < A < 1 

(4.6) \\dxd,Axm^,<CRQ^/^{t). 
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Now, we observe that by Jensen's inequality, 



/ \d,Ai{t,x)~d,A2{t,x)\ 

JM3 



dx 



dxd.,Ax{t,x)d\ 







dx 



< / / \dxd,Axit,x)\UxdX 



< sup / \dxdxAx{t,x)\'^dx, 

0<A<1 JR3 



and similarly for ||Ai(i) — ^2(^)11^2. Then, we use (|4.5p and ()4.6p to derive the 
estimate 

(4.7) \\A,{t) - A2(i)||^.(s^) + ||9,Ai(i) - d,A,{t)\\^, < CflQi/2(t). 
Finally, we combine pTTj) . (|i?^ and (|i?71) to conclude that 

« < C,0(,). 

Uniqueness then follows as a trivial consequence of Gronwall's lemma. D 

Proof of TheoremUl Let /o G Co'"(M^_K), /o > 0. Fix Xq > and Pq > such 
that fo{x,p) = for \x\ > Xq or |p| > Pp. We introduce the following iterative 
scheme. For i e / and z — {x,p) € R'^ x R-^, define 

/0(i,z)-/o(z). 

For n G N, assume that /" : / x R® ^- R is given and define 

(4.8) ^-it,x) ^ f f f^t,y,p)^^ 

JR3 Jm^ \y - x\ 

(4.9) A^it,x) = \[ f [±d + Lo(g>Lo]vAj''{t,y,p)-,'^^'^y 



2 Jr3 Jr3 ly - x\ 

where 



l+lp-Anf 



Denote by Z„ = (X„,P„)(s, t, z) the solution of the characteristic system 

(4.10) X„{s,t,z) = WA„(s,X„(s,i,z),P„(s,t,z)) 

(4.11) P„(s,t,z) = -[V$"-t;^„VAy (s,X„(s,i,z),P„(s,t,z)) 

with Zn{t, t, z) — z. We define the {n + l)-th iterate of the distribution function by 

r+i(i,z) = /o(Z„(0,i,z)). 
For convenience we shall also define the sequences 



p^(t,x)= i f\t,x,p)dp and i„(i,a;)=/ VA„f"'it,x,p)dp. 

Step 1: In view of the Lemmas and Remarks in Section [2l and Lemma [7] in Section 
[3 the sequence {(/",$",A„)} is wch defined. In particular, /" G C^{I,C^'°'(R%R), 
/" > 0, and ($", A„) G (7^(1, C2'"(R3); R x R^). For each n, the regularity in time 
of the potentials is the one of /". This is trivial for $". As for A", see Remark [2] 



llp"WllLo.<-^||/olLson?(i). 
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For t e I set Po{t) = Pq and for n e N define 

P„(i) = {\p\:30<s<t,xGW':r{s,x,p)^0} 
= {|F„_i(s,0,z)| : 0< s < t,ze supp/o} 
It is clear tliat supp/"(i) C {(x,p) e R^ ^ ]r3 . |^| < Xa + t, \p\ < P„(i)}- Also, 

\\rmLi = \\M\Lv i<9<oo' *e/, nen 

and we have the estimate 

3' 
Since \jn\ < |/5"|, the known estimates on the potentials imply that 

\\'^''it)h^ +\\A^{t)h^ <C{fo)Pn{t), 

and finally 

(4.12) ||9,$"(i)|j^^ + \\d,A,M\L^ < C{h)Pl{t). 

Step 2: For some T > there is a non-negative, non-decreasing V G C([0,r[;M) 
depending on the Cauchy datum only, such that for all n e N U {0} and < t < T 

Pn{t) <V{t). 

Indeed, for n e N the characteristic equation (|4.11l) and the estimate (I4.12p imply 

(4.13) < Po + C(/o) / Pl{T)dT. 

Let T > be the life span of the solution of the integral equation 

(4.14) V{t)=Po + C{fo) f rHT)dT. 

Ja 

Hence, Po(i) < V{t). Suppose P„(t) < V{t) for some n G N. Then, in view of 
(|4.13l) . this estimate also holds for P„+i(i), which proves the claim. As a result, all 
estimates in Step 1 are uniform in n on any subinterval [0, T] C [0, T[. In particular, 
for all n e N and < i < T, we have 

(4.15) ||p"(t)|L. + ||jn(i)|Lo. + \\d.^'\t)\\^^ + Wd^AniML^ <C^^ C(t,/o). 

X X X X -^ 

For future use, we notice that the maximal solution of (|4.14p is given by 

(4.16) V{t)^Po{l-C{fo)Poty\ 0<t<T={C{fo)Po)~\ 

withC(/o)=3(27r)2/3||/o||V3 Wf^f/^ , 

Step 3: We claim that for every fixed < T < T 

(4.17) ||a.p"(i)iL. + ||a.j„(i)||^^ + \\dl'^-{t)\\^^ + \\dlA,,{t)\\^^ < c^, 

for all n e N and < i < t. 

To start with, we estimate the space derivatives of the characteristic curves. To 
ease notation, we write (X„,P„)(s) = (X„,P„)(s,i,a;,p). Recall 

VaJs,X„{s),P„{s)) = v{Pn{s),Anis,X„{s))), 
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where v is C^ in its argument. Hence, since (by abuse of notation) we have 
{dccXn{t),dxPn{t)) = (1,0), the uniform bounds in (|4J5l) lead to 

\d,Xjs)\ < \d^Xjt)\+J |a,[t;(P„(T),A„(T,X„(T)))]|dT 

^0 



< 1 + C^ / i\d,X„{r)\ + \d,Pn{T)\) dr. 
Jo 



Similarly, 



\d.Pn{s)\ < \d,Pnit)\+ f \d^[V^-v\yAl^]{T,X^a{T),Pn{T))\dT 

J s 

< C° /* (l + ||a^$"(r)||^^ + \\dlAUT)\\^^ 
Jo "" "" 

X ( |a,x„(T)| + |9,p„(T)| 



dr. 

These two estimates and the Gronwall's lemma yield 

|a,X„(s)| + |9,F„(s)| 

< 



exp|c°^ (l + ||a^*"(r)|L^^ + ||5|A„(r)||^^)dr 



As a result, we also have 

|5.p"+i(i,x)| < / |9,[/o(^„(0,t,a;,p))]|dp 



bl<-P(t) 

t 



Similarly, after using the product rule and the known estimates, there exists a 
sufficiently large constant C^ such that 

\dxjn+lit,x)\ < / \dx[v{p,An+l{t,x))fo{ZniO,t,X,p))]\dp 

J\p\<V{t) 

< C^expjc^^ {i+\\dl<^"ir)\\L^ + \\dlA^ir)\\^^)dTY 
Hence, in view of Lemma [6l we have for all < i < T that 

< C^{l + ln+ \\dxp"+\t)\\^^ + ln+ \\dxjn+i{t)h^) 

< C7?+C5 r(||a2$"(r)||^^ + ||92A„(r)||^^)dr. 

Jo "" "" 

Since the right-hand side is bounded for n = 0, induction in n yields 
\\dl^'\t)\\^^ + ||9X(i)|Lo. < C^expjC^T} , 

for all n S N and < t < T. In turn, this provides a uniform bound on the 
derivatives of the iterates for the current and density functions. 
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Step 4: We show that {/"} is Cauchy in the uniform norm on [0, T] x M^. To start 
with, notice that 

\r+\t,z)-r\t,z)\ = i/o(z„(o,t,z))-/o(z„_i(o,i,z))i 

(4.18) < C|Z„(0,i,z)-Z„_i(0,i,z)|. 

On the other hand, by using the estimates in the previous steps, it is not difficuh 
to check that the characteristics equations lead to 

|X„(s) - X„_i(s)| 

< f |l>(P„(T),A„(T,X„(T)))-«(F„_l(T),A„_l(T,X„_l(T)))|dr 
J s 

< C I (|X„(t) - X„_i(t)| + |P„(t) - P„_i(t)| 

+ ||A„(r)-A„_i(T)||^<^)dT, 
and 

|P„(s) - Pn-l{s)\ 

< J ( |V$"(r,X„(r))-V$"-i(r,X„_i(r))| 
+ 1 (t;^„VA;)(r,X„(r),P„(r)) 

- («A„_,VA^_i) (T,X„_i(r),P„_i(r))|) dr 

< C /" ( |X„(t) - X„_i(r)| + |P„(r) - Pn-iiT)\ 

a; 

+ ||A„(r) - A„_i(r)||^^ + ||a:,A„(T) - a:,yl„_i(r)||^^) dr. 

X X J 

Therefore, after adding the above expressions, Gronwall's inequahty yields 

|z„(o,i,z) - z„_i(o,i,z)|<c/ (||a,$"(T)-a,<i>"-i(T)||^„ 

Jo ^ 

(4.19) + ||A„(T)-A„_i(T)||^^^ + ||9,A„(T)-a,A„_i(T)||^^^)dr. 

Now, to produce a Gronwall's inequality resulting from (I4.18P and (j4.19p . we 
look for suitable estimates on the right-hand side of (|4.19p . To start with, let 
R = max {Xo + T, V{f)]. For all n e N and < t < T we have 

supp/"(t) C Bh X Br. 

Linearity and Lemma [5] yield 

\\d.<^'\T) - d,^-\T)\\^^ < C ||p"(r) - p"-^(r)||^f ||p"(r) - p''-\r)f/^ 

(4.20) < Ci, II r(r)-r-i(r) 11^^^^. 

To estimate the terms involving the vector potential, we proceed as follows. Ac- 
cording to the definition of the iterates, it is clear that for each n G N they satisfy 

dtr+^ + VA^ ■ v,r+i - [v$" - v\yAi;\ ■ Vpr+^ = o. 
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!A„(r) - A„_i(r)||^.(^^) < Cr ||r+i(r) - /"(r)] 



Hence, by the uniqueness proof now in terms of the iterates, (see (14. 7p ). 

which can in turn be estimated in terms of 11 /"^^(t) — /"(''■) 11 roo • 

Also, if we write respectively /" — /"^^ and An — An-i in Lemma [9] instead of 
dtf and dtA, it is easy to check that 

P„(r) ~ An-,{T)\\^^ 

Similarly, after a slight modification of the proof of Lemma[TO] (it is actually simpler 
here since the Vlasov equation is not used), we find 

ll5,A„(r) - a,A„_i(T)||^^ 

a; 

< Cn (||r(r) - r-\r)\\^^^ + ||(A„(r) - A„_i(r))||^.(^^)) . 

Therefore, the previous three estimates yield 

||A„(t) - A„_i(r)||^^„ + ||a,A„(r)-a,A„_i(r)||^^^ 

(4.21) < Cr (||r+\i) - r 11^^ + ||r(T) - r-\T)\\^^ ) . 

Hence, if we combine (|4.18p and (14.19^ with (|4.20l) and (|4.2ip . a use of Gronwall's 
lemma gives 

ll,r+i(i)-r(t)ll^^ <c [ \\r{r)-r-\T)\\^^ dr, 

which by induction, readily implies the claim. It follows that {/"} converges uni- 
formly to some / e C([0, T] x R^. ]r) ^nd for all < i < T we have 

supp/(t) C Br X Si^. 

Finally, if we respectively define p, $ and A according to ()1.4p , p.2p and ()3.3p , we 
have that p, $ and A are Cb, and p" — >• p, $" — > $ and A„ — > A hold uniformly 
on [O,?] X R^. The latter follows from (|4.2ip . The uniform limits va^ — > va and 
'^A„f" -^ ua/ can be easily checked, and therefore j„ — > JA uniformly on [0,T']xR'^, 
with JA e Cb{[0,f] X R3;R3) defined by (fTAll . 



Step 5: Actually / G C^(/xR^;R), as we shownext. Indeed, in view of Step 4 and, 
respectively, (I4.20p and (I4.2ip . the sequences {dx^"} and {SjjAn} are uniformly 
Cauchy on [0,T] xM.^. Moreover, by Lemma [6] we have 

\\dlA„,{t) - dlAr,{t)\\ < C \r-' ||j„,(i) - jn{t)h. 

+ h\\dxjm{t) - dxJn{t)\\L^ + (l+\n{R/h) \\jrn{t)~Jn{t)\\L^) 

and similarly for 9^$"(t). Hence, the known estimates and the fact that we can 
choose h arbitrary small imply that {9^$"} and {9^A„} are uniformly Cauchy on 
[0, T] X R'^ as well. Therefore, we have (with a slight abuse of notation) 

($,A), a,($,A), a2($,A) eC([0,t]xM-^RxM-''), 
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and so the characteristic flow Z e C^([0,T] x [0,T] x R^) induced by the hmiting 
field is in turn the hmit of the sequence {Zn}- As a resuh, the function 

/(i,z)= hm /o(^„(0,i,z)) = /o(Z(0,i,z)) 
has the claimed regularity and the triplet (/, $, A) satisfies p.l|) - p.4|) . 



Step 6: We show that the potentials $ and A have the required regularity in time. 
Define fx = A/" + (1 - A) /", < A < 1. This definition is just like the one in the 
uniqueness proof but in terms of any two elements of the sequence {/"}. Let A\ 
be the vector potential induced by f\. Following the lines in the proof of Lemma 
|S1 it is not difficult to check that (this is analogous to (I3.1ip ) 

\d\dtdxAx\^ = fxdxdtAx-dxdtVA^ 



Br J Br 

dxdtA ■ [dxVA^dtfx + dtVA^dxfx + VA^dxdtfx] 

Br J Br 

(4.22) + dxdt f f i (V • Aa) (V • dxdtjx) ■ 

Jk^ J Br r 

Since V • Ax ~ 0, the third integral in the right-hand side vanishes. On the other 
hand, by using the notation of the Appendix, we have Oxvax = —DvA^dxAx, also 
dtVAx = -DvA^dtAx, and 



dxdtVA, = -DvA.dxdtA - D^VA.dxAxdtA 



X- 



Therefore, since by Step 5 \dtfx\ < Cr, and by Corollary [2] I^jAaI < Cr, we 
obtain from (|422|) that 

\dxdtd,,Axf + f [ -JL= (\dxdtAxf ~ \VA, ■ dxdtAxf) 

J Br J Br y/l+gi^ J 

< Ch f f \dxdtAx\[\dxAx\ + \dxfx\ + \dxdtfx\] 

J Br '^ Br 

(4.23)< CR\\dxdtAx{t)\\^, \\\dxAx(t)\\^,^SR) + \\dxhmLs^ +\\dxdtfx{t)\ 



i; 



In the last step we have used the Cauchy-Schwarz inequality. Now, define 



G^ 



= sup (\\nt)-rmL^ +\\dtnt)-dtrmL^) 



mn 

0<t<T 






which in view of Steps 4 and 5 converges to zero as n, to — >■ oo. If we use the estimate 
(|4.5p for the iterates, i.e. \\dxAx{t)\\^2/B^ < CR\\dxfx{t)\\i^2 , we find that the 
expression in square brackets in the right-hand side of (|4.23p can be estimated as 

\\dxMt)\\Ll(BR^ + \\dxfx{t)h^.^ + \\dxdjx{t)h 

< Cr (\\dxfx{t)h2 +\\dxfx{t)h^ +\\dxdtfx{t)h^) < CRG,^n, 
\ x,p s,p ^,p/ 

uniformly in A. On the other hand, since \vA)^ \ < 1 strictly, we can reason as in the 
proof of Lemma [8] to find a lower bound on the left-hand side of (|4.23p . This lower 
bound can then be estimated as 

Wdxdtd^Aximli + \\pTit)dxdtAx{t)f^^^^^^ < Cr WdxdtAxmLiiBR) G„™. 
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Consider the first terra on the left-hand side. Poincare's inequality and the above 
estimate imply that 

Then, the last two estimates yield 

(4.24) \\d^dtAx{t)\\^,^s^^ + Wdxdtd^A^m^, < CRGmn- 
On the other hand, by the definition of A\, we have 

dxdtAx{t,x) = j j lC{x,y)dxdt[vAxf\{t,x,p)]dpdy. 

^ Br J Br 

Therefore, after taking the product rule in the integrand, we may proceed as in 
Lemma ini to obtain the estimate 

(4.25) WdxdtAxim^^ < Cr (g„„, + \\dxdtAx{t)\\^,^^s^^ , 
where again we have used \dtfx\ < Cr and 19*^^1 < Cr. Similarly, since 



dxdtdxA\{t,x) ^ I I dx)C{x,y)dxdt[vAxfxit,x,p)]dpdy, 

J Br ^ Br 

we can proceed as in Lemma [TU] to find 

(4.26) WdxdtdxAxmL^ < Cr (G™„ + \\dxdtAx{t)\\^^^g^^) . 

Hence, since \dtAjn — dtAn\ < /q \dxdtAx\dX < aupx\dxdtAx\ and similarly for 
\dtdxAjn — dtdxAjil, we can gather the above estimates to find that 

WdtA^it) - dtA^m^^ + \\dtdxA^{t) - dtdxA^{t)\\^^ < CRG^n. 

Therefore, the sequences {9f A„} and {dtdxAn} are uniformly Cauchy and we have 
that dtAn — >■ dtA and dtdxAn — J> dtdxA uniformly on [0, T] xR^. In turn, the former 
limit and Steps 4 and 5 imply the uniform convergence dt{vA„f^'') — >■ dtivA.f), and 
so dtjn — >■ dtJA- Also, 9tp" — )■ dtp. Hence, just as in Step 5, the sequences {9i9^$} 
and {dtd^A^ are uniformly Cauchy on [0,T] x M."^. Therefore, since trivially 9t$ 
and dtdx^ are continuous on [0,T] x R^, we conclude that 

dt{<i>,A), dtdx{<^,A), dtdl{<i>,A) e C{[0,f] X r3). 

Having proved the claim, and since < T < T was arbitrary, we conclude that 
/ G C^([0, T[xR^; R) is a classical solution of the relativistic Vlasov-Darwin system. 

Step 7: Moreover, f{t) G C^^"(R6;R), < a < 1, for each < t < T. In 
view of Remarkm this holds if ($, A) (i) e C2^"(R3; R x R^). But, since we have 
iP,.u)(.t) e C^{M.^; R X R3) c C^(R^; R x R^), the regularity needed for the poten- 
tials is guaranteed (see the last lines in the proof of Lemma [7]). 

Step 8: The proof of the continuation criterion is as follows. Let / be the solution 
of the RVD system previously obtained, which clearly satisfies /Ij^q' ^^ shown in 
(I4J61) . the life span of / is T = (C(/o)Po)"^ with 

C(/o)- 3(2^)^/3 ll/oll^/^ ll/olli/i . 

x,p x,p 

Define Pt = sup {\p\ : 30 < t < T, x £ R^ : f{t, x,p) ^ O} and assume that Pt < oo 
but T < oo. We claim that this is a contradiction. 
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Fix < to <T and consider /(to) as a Cauchy datum of the RVD system, which 
is guaranteed by Step 7. Known estimates yield 

ll/(io)ILi = ll/olLi , ll/(to)ILoo - ll/olLo. . 

j;.p .r.p x,p a:,p 

Thus, C{f{to)) — C(/o). Define e = (C(/o)Pt) , which does not depend on to- 
Steps 1-3 imply that all uniform estimates on the sequence of approximate solutions 
induced by /(to) hold on [to, to + e[. Then, /(to) yields a unique classical solution 
of the RVD system on that interval. 

But we could have fixed to arbitrary close to the life span T < oo of / and so 
extend this solution beyond T, which is a contradiction. Hence, we have shown 
that Pt < oo implies T — oo. This, and the uniqueness result, conclude the proof 
of Theorem m D 

4.2. Global Solutions. If additional conditions are imposed on the Cauchy datum 
in Theorem [U then the local solution found in the previous section can be extended 
globally in time. We prove this result next. We start by defining the set where the 
Cauchy datum will be taken from. For Xq > and Po > given let 

V = {f eC^-°'{R^;R), 0<a<l: 



/ > 0, WfWwt- ^ 1' ^^PP/ C ^Xo X Bp,] 



Theorem 2. There exists a S > such that, if fo (zV with ||/o||roo < S, then the 
classical solution of the RVD system i3.1\) - (^^^ with Cauchy datum /o is global in 
time. Moreover, for t > this solution satisfies the decay estimates 

(4.27) \\pmLr + \\jA{t)h^ < Ct-^ 

(4.28) \\d.Ht)\\^^ + \\d,AmLr < Ct-' 

(4.29) \\dmt)\\^^ + \\dlAit)\\^^ < Ct-Mn(l + t) 



We first introduce some technical results and postpone the actual proof of The- 
orem [2] to the end of this section. The following lemma shows that a sufficiently 
small Cauchy datum leads to a classical solution of the RVD system which exists 
on any given time interval and induces potentials whose derivatives can be made 
as small as desired. 

Lemma 11. Fix e > and T > 0. There exists 5 = <5(e, T) > such that, if fo G T> 
with ll/oll roo < <^; then the classical solution of the RVD system with Cauchy datum 

a:,p 

/o exists on the time interval [0,T] and induces potentials satisfying 

(4.30) \\dtA{t)\\^^ + \\d,A{t)\\^^ + ||a,$(t)||^^^ + \\dlA{t)\\^^ + ||9^'i>(t)||^^^ < e 

for allQ<t<T. 

Proof. In view of Lemma [6l and since \ja\ < p and 

\\dxPit)h^ + ||9.M(t)ILoo < C^, < t < T 

hold (the latter proved just as in Step 3 in Theorem [Jl with the estimates applied 
to the solution instead of the iterates) , the space derivatives of A satisfy the same 
estimates as the space derivatives of $. Hence, the proof is mutatis mutandis 
the proof of [HI Lemma 4.2] for the Vlasov-Poisson system, as far as the space 
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derivatives of the potentials are concerned. As for dtA, the result follows suit in 
view of the estimates in Lemmas |8] and 13 D 

To proceed, we now define the so-called free streaming condition for classical 
solutions of the RVD system. 

Definition 4. Fix /3 > and a > 0. A classical solution of the RVD system is said 
to satisfy the free streaming condition of parameter (3 (FSj3) on the time interval 
[0, a], if it exists on [0,a] and induces potentials satisfying the estimates 



\\dtA{t)\\^^ + \\d.A(t)\\^^ + ||9.<i>(t)||^^^ 


< /3(l+t)-3/2 


dlA{t) ^^ + dl^{t) ^^ 


< /3(f + t)-5/2 


for alio <t<a. 





Lemma 12. There exists S > 0, f3 > and a positive C — C{Xq,Pq) such that 
any classical solution f of the RVD system having a Cauchy datum foCzV with 
ll/oll Too < (5 and satisfying (FS/3) on some interval [0, a], also satisfies the estimates 

(4.31) \\d^<^{t)\\L^ + \\d.A{t)\\^^ < Ct-\ 

(4.32) \\dlHt)\\^^ + \\dlA{t)\\^^ < Ct-'^lnil + t), 
for all < t < a. 

Proof. By virtue of Lemma O the following estimates hold 
\\dMt)\\^^ + \\d,Ait)\\^^. < CWfoC \\pit)fi!', 
\\dl^{t)\\^^ + \\dlA{t)\\^^ < Ct-^\\\fo\\]!' +\\dMmL^ + \\d.jA{t)h^ 






{l + \nt*)t^\\p{t)\\ 



L° 



t > 1, 



where the latter is a consequence of setting R = t and h = t^^ < R in the cited 
lemma. We claim that for some suitable constant C = C{Xo,Po) > the charge 
and current densities satisfy 

(4.33) \\pmL^ + \\jA{t)h^ < Ct-\ 

(4.34) \\d.p{t)\\L^ + \\d,umL^ < C. 

If true, then the lemma follows. To prove the claim, we first introduce some tech- 
nical results which we present as a sequence of steps. 

Step 1: Let < s < i < a. Denote by {X,P){s) = {X,P){s,t,x,p) the solution of 
the characteristic system 

X{s) = va{s,X{s),P{s)) 

P{s) = -[^^ + v\VA^]{s,X{s),P{s)), 

with {X,P){t) = {x,p). Denote also Dva{s) ^ Dva{P{s),A{s,X{s))), where the 
matrix Dva is as given in the Appendix. Consider the system 

^(s) = dpXis)^{s-t)DvA{t) 
r]{s) = DvA{s)dpP{s)-DvA{t). 
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Notice that ^{t) = r){t) = 0. We show that for some C = C{Xo,Po) > 

(4.35) \as)\<fiCe^^{t-s). 

Indeed, on the characteristic curves, we have 

iis) = dpX{s)-DvA{t) 

- DvA{s)[dpP{s) - d^A{s,X{s))dpX{s)]- DvA{t) 
= 7j{s) - DvA{s)d^A{s, X{s)) [e(s) + {s- t)DvA{t)] . 

Therefore, since |_Di;/i(s)| < C, a use of (FS/3) yields 

leWI < f\viT)\dr + (3cf\l + r)-'/'mr)\ + it-T)]dr 

J S J S 

(4.36) < pCeP''l{t-s)+ f\vir)\dT), 



where the Gronwall's incquahty has been used in the last step. 
On the other hand, we have 

r,{s) = DvA{s)dpP{s) + D\a{s) [p{s) - A{s, X{s))\ dpP{s) 

In view of the characteristic system, it is not difficult to check that 



+C\\d,A{s)\\^^\dpP{s)\. 



Hence, since A* — dgA^ + va ■ VA*, i — 1,2,3 and |D^u^(s)| < C, the above 
inequality and (FS/3) yield 

\ri{^)\ < [\\dMs)\\^^ + \\d',A{s)\\^^+\\d^A{s)\\l^)\dpX{s)\ 
+ (l|9.$(.s)||ijo + \\d^A{s)\\^^ + \\dsA{s)\\^^) \dpP{s)\ 
< 2/3 (1 + s)-"/^ \dpX{s)\ + /3 (1 + s)-'/^ \dpP{s)\ . 

Now, by the definition of ^(s) and ri{s), we have |9pX(s)| < \(,{s)\ + C {t — s) and 
\dpP{s)\ < C (|??(s)| + 1), the latter as a result of |Dw^"'^(s)| < C, as it can be easily 
checked. Then, Gronwall's inequality implies 



\rj{s)\ < (3Ce 



I3C 



{l + ry'/'\^{r)\dr 



(4.37) 



{1 + T)-''^\t - t) + {1 + t)-^^^' 



dr 
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Both (|436l) and (j437l) then lead to 



|?(s)| < /SCe'^M (t - s) + 



fj\l + a)-"^\a<y)\dadr 



(1 + gY^''^ {t-G) + {l + (jy^'^] dadA 






il + a)-^^^\^ia)\dTd<J 



S -/ S 

PC 



■t 



-5/2 (-+_ ^) + (1 _L^\-3/2 
-3/2 

-3/2 /-+ _ „\ I n ^ ^^-1/2 



drdcr I 



< pCe^^ [{t-s)+ I {l + ay'Ua)\dcj 



{l + a)-''''{t-a) + {l + ay 



da 



Finally, since the last integral is less than 3 (t — s), another use of Gronwall's in- 
equality yields (|4.35p . 

Step 2: For /3 > small enough, there exists a C = C{Xo, Pq) > such that the 
mapping ^^^(0, t, x, •) : R"^ — ?► R"^ has Jacobian determinant satisfying 

\detdpX{0,t,x,p)\>Ct^, < t < a, x e M^ peR^. 

For t = this is obvious. Let < i < a. Without loss of generality, we shall assume 
that < /? < 1/2. Then, by the characteristics and (FS^) we have 

_ ft 

(4.38) \Pit)\<Po + l3 (l + T)"'/'dr<Po + l. 

Jo 

Also, in view of the estimate on the vector potential given in LemmalBl and recalling 
that /o € T), is not difficult to check that 

Denote g = \p — A\. Hence g < C{Xo,Po) and therefore the relativistic velocity 
satisfies |z;^| < t^ < 1, where ly depends only on Xq and Po. Now, we have that 

DvA = (1 + g^y^^^ [id -VA(^ va]- Then, since by Step 1, |C(0)| < l3Ce'^^t with 
^(0) = dpX{Q) + tDvA{t), we have for some /3 > small enough that 



\A+? 



t 



dpX{0) + id 



yi+? 



t 



e(0) + va(E)va 



(4.39) < /3Ce^'^ + iy = 7 < 1. 

Therefore, a positive constant C = C{Xo, Pq) exists such that 



\detdpXiO)\ 



(1+5^) 



2^,3/2 



det 



t 



dpX{0) + ±d-±d 



>Cf^ 



Step 3: For every < t < a and x e M^, the mapping X{0,t,x, •) : R^ -> M^ is 
bijective. Indeed, for p,q £ R^, let 

Px = Xp+{l-\)q, gx^g{t,x,px)^\px-A{t,x)\, < A < 1. 
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In view of (|4.39p in Step 2, we have 
\XiO,t,x,p) - XiO,t,x,q)\ 

r-l 



dpX{0,t,x,px){p-q)dX 



-id+ id 



ViT^ 



t 



dpX{0,t,x,px) 



t {p ~ q) 
v/1 + 5!^ 



dX 



, [^ dX , , /"^ dX 



/o y 
> (l~7)b-(7|i 



which shows that the mapping is injective. It is also surjective, since the open range 
X{0, t, X, M.^) —M.^. If not, there exists a boundary point xq so that X{0, t, x,pn) — > 
xq ^ X{0, t,x,R'^) as n — > 00, for somep„ -> po G ^^- By continuity X{0, t,x,po) = 
Xq, which is a contradiction, and the assertion fohows. 



Step 4: Then, Steps 2 and 3 imply that the mapping X{0,t,x,-) : M'' ^- M"^ is 
a C^-diffeomorphism. In particular. Step 2 implies that for some constant C = 
C{Xo, Po) > 0, the inverse mapping ^-^(O, t, x, ■) -.W^ ^ R^ defined by X k^ p{X) 
has Jacobian determinant satisfying 



\detdpX-\0,t,x,p{X))\ < Ct-^, < t < 



a, x e 



We can now deduce the estimates (|4.33[) - (|4.34p for the charge and current densities. 
Indeed, bearing in mind that /o G D, we have 



p[t,x) 



foiXiO,t,x,p),PiO,t,x,p))dp 
fo{X,P{0,t,x,p{X)))\detdpX-\0,t,x,p{X))\dX 



< Ct 



-3 



where C = C(Xo,Po) > 0. Then, since \ja\ < p, (|i35)) indeed holds. 

To prove (|4.34p we proceed as follows. In view of (FS/3) with ji = 1/2, and 
recalling that |92;Wa| < C \dxA\ and /o G T), we have that 



(4.40) ||a,p(t)||^ 

(4.41) 

and 

(4.42) 



< C{Po + iy\\dJ{t)\\^^ 

< C(Po + l)'(||9.A(t)||^^^||/o||^ 

< c(Po + i)'(i + l|a,/(i)||i. 



\d.fm 



\dxfit,x,p)\<\d.xPiO,t,x,p)\ + \dxX{0,t,x,p)\. 



Hence, the proof will be completed if we provide a uniform bound on the space 
derivatives of the characteristic curves. Similar to the computations in Step 1, it is 
not difficult to check that for < s < t < a 



\dxX{s)\ <1 + C f {\dxP{T)\ + \\d.,A{T)\\^^ |9.X(r)| 



dr, 
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and by (FS/3) and Gronwall's lemma 

(4.43) \d^X{s)\<c(l + J \d^P{T)\dT 

Also, 

\d.P{s)\ < C [ {}d'M')\\L^ + PIMs)\\l^ + \\d^.A{s)\\l^) \d^X{T)\ dr 

+ C j \\d^A{T)\\^^\d^P{T)\dT 

(4.44) < C j {l + T)-'''^\d^X{T)\dT. 
Therefore, both (|i^ and (|TiI)) yield 

\d^X{s)\ < C + C I f {l + ay^/^\d^X{a)\d(jdT 

< C + C f f {l + ay^^^\d^X{<7)\dTda 

< C + C f {l + ay^^^\d^X{a)\da. 

Gronwall's lemma then provides a uniform bound on |92:-'^('S)|, which in turn pro- 
duces a uniform bound on \dxP{s)\ via (I4.44p . As a consequence 

\dxX{0,t,x,p)\ + \dxP{0,t,x,p)\ <C, < t < a, a;eR3, pe]R^ 

which implies (|4.34l) via (|4.40l) - (|4.42l) . This concludes the proof of the lemma. D 

Proof of Theorem \^ By virtue of the previous lemmas, the proof is almost 
identical to the proof of [131 Theorem 4.1] for the Vlasov-Poisson system. 

Indeed, let /3, 5 > and C = C(Xo,-Po) > be suitable for Lemma [T^] to hold. 
Fix To > 1 such that for all t > Tq 

(4.45) Ct-^ < - (l + ty^^\ C7(l + lni)i-3 < ^(l4_t)-5/2^ 

Now, by letting 5 > be smaller if necessary. Lemma [TT] implies that the Cauchy 
datum /o G 2? with ||/o|lioo < S yields a classical solution / of the RVD system on 
the maximal existence interval [0,T[ with T > Tq, and 

\\dtA{t)\\^^ + \\d,A{t)\\^^ + ||a.ci>(t)||^^^ 

+ \\dlA{t)\\^^ + \\dlHt)\\^^ < ^ (1 + To)-^/^ , 

for all < i < Tq. Hence, / satisfies the free streaming condition (FS/3) on [0,To]. 
In fact, the continuity of the left-hand side of the above inequality implies that there 
exists a maximal Tq < Ti < T such that / satisfies (FS/3) on [0,Ti[. Therefore, 
Lemma inland ^fAM\i imply that for all To <t<Ti 

\\dx^t)\\^^ + \\dxA{t)\\^^ < ct-2 < ^(i + i)-3/2^ 

\\dl^[i)\\^^ + \\dlA{t)\\^^ < C(l + lnt)t-3 < ^(l + t)-5/2. 
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Then, a continuation argument yields Ti = T, and by (j4.38|) . we deduce 

sup{|p| ■.30<t<T,xeR^ : f{t,x,p)^0} <Po + l. 

Therefore the continuation criterion in Theorem [T] imphes that T — oo, and thus 
the solution / is global in time. The proof of Theorem [2] is complete. D 

Appendix 

T. P~ A id-VA(E>VA ^ 
boi va = — , set -UvA = — — OpVA- 

^l + \p-A\' ^l + \p-A\' 

Then O^va — —DvAd^A and dfVA — —DvAdfA. Clearly, |I?w^| < C, and so 
\d,VA\ < C\d,A\ and \dtVA\ < C \dtA\. Also, \d^VA\ < C; \d,dpVA\ < C \d,A\; 

\dlvA\ < c{\dlA\ + \d,A\'y, \d!vA\ < c{\dfA\ + \dtA\^y,\dtdpVA\ < C\dtA\ 

and finally \dtd^VA\ < C i\dtd,A\ + \dtA\ \d,A\). 
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